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ABSTRACT 



Vh ' We study numerically the magnetic susceptibility of the hierarchical model 

c^ ' 

with Ising spins (a = ±1) above the critical temperature and for two values of the 

epsilon parameter. The integrations are performed exactly using recursive methods 

which exploit the symmetries of the model. Lattices with up to 2^^ sites have been 

used. Surprisingly, the numerical data can be fitted very well with a simple power 

law of the form {l—(3/(3c)^'^ for the tt;/io/e temperature range. The numerical values 

for 7 agree within a few percent with the values calculated with a high-temperature 

expansion but show significant discrepancies with the epsilon-expansion. 



1. Introduction 

The renormalization group (RG) method is a powerful tool to handle criti- 
cal phenomena and to approach the continuum limit of lattice models. However, 
its practical implementation usually requires approximations. In his original pa- 
per, Wilson made order of magnitude estimates of various terms contributing to 
the partition function of the Landau- Ginzburg model and derived the so called 
approximate recursion formula. In this approximation, the RG transformation is 
reduced to a single integral equation which can be studied using numerical methods 
or functional analysis. 

Recursion formulas closely related to the approximate recursion formula hold 
exactly for the hierarchical models. This is due to the large group of symmetries 
of the hamiltonians of these models. The RG transformation for these models has 
been studied in great detail and rigorous results concerning the epsilon expansion 
of the critical exponents are available in the literature. 

The fact that the RG transformation can be handled easily for hierarchical 
models suggests the use of these models as an approximation for nearest neighbor 
models. The main technical problems in proceeding this way are: how to derive 
explicitly the approximate models and how to improve systematically the approxi- 
mation. Recently, one of us '' has answered these questions for the gaussian models 
where everything can be calculated explicitly. In order to extend this method to 
interacting models, one should be able to calculate the average value of perturba- 
tion terms added to the hierarchical hamiltonian. For Ising models, where the spin 
a takes only the values ±1, this task can be carried numerically in an efficient way. 
It has been suggested that such an approach might shed some light on Polyakov's 
conjecture for the 3D Ising model. In preparation for these calculations, we first 
checked the agreement between the numerical calculations for the (unperturbed) 
hierarchical Ising model and analytical results. In doing so, we found surprising 
results which are reported in the following. 

In this paper, we calculate numerically the magnetic susceptibility per site - the 



susceptibility for short - of the hierarchical Ising model as a function of the tem- 
perature and for two values (0 and 1) of e, the parameter used in the e-expansion. 
Calculations have been carried with up to 2 sites. The numerical integration 
made use of the symmetries of the model in order to cut down the time of com- 
putation logarithmically. However, no approximations have been made and the 
numbers shown below are exact up to round-off errors. These errors were ana- 
lyzed by changing from simple to double precision. In all the cases considered, this 
only affected the fifth significant digit of the susceptibility in the worse cases. Our 
calculations have been mostly restricted to the high temperature region and its 
boundary. In other words, the parameter j3, proportional to the inverse tempera- 
ture, will run between and a critical value (3c- However, at the beginning, a few 
calculations will be made in the low temperature region in order to locate f3(,- 

Surprisingly, we found that the numerical data can be fitted very precisely 
with a simple power law of the form (1 — (3/f3c)~^ in the whole high-temperature 
region. As a consequence, it is possible to calculate with good accuracy the free 
entries of this parametrization, 7 and j3c, using the two first coefficients of the 
high-temperature expansion of the susceptibility. On the other hand, the values 
of 7 obtained numerically, differ significantly from the values obtained in the e- 
expansion . This does not mean that the approximations made to calculate the 
e-expansion are incorrect, but only that for the models considered here, we can 
never get near the region (in the space of theories) where these approximations 
can be made. 

This paper is organized as follows. In section 2, we introduce the hierarchical 
model and the numerical method of integration. In section 3, we check the scaling 
laws for the variance of the total spin and we determine (3c- In section 4, we analyze 
the temperature dependence of the susceptibility and we show that the data can 
be fitted very well with a simple power law. In section 5, we compare the results 
with the high temperature expansion and the e-expansion. Finally, we discuss our 
present understanding of the results in the conclusions. 



2. The Hierarchical Ising Model 

In this section we describe the hierarchical Ising model and the basic ideas of 
the numerical calculation performed. Hierarchical models are specified by a non- 
local hamiltonian bilinear in the spin variables and a local measure of integration. 
We consider here the case of a Ising measure, where the spin take only the values 
±1. The hamiltonian of a hierarchical model with 2"" sites can be written as 



1 " 

^E(f)' E ( E -(.. ..))' (2-1) 



H - 

2^^M 

I— I in,...,ii+i ii,....,ii 

For convenience, we have labeled the sites with n indices in ii, each index being 

1 or 2. In order to visualize the meaning of this notation, one can divide the 2" 
sites into two boxes, each containing 2"""^ sites, li in = I, the site is the first box, 
if in = 2, the site is in the second box. Repeating this procedure n times (for 
the two boxes, their respective two sub-boxes, etc.), we obtain an unambiguous 
labeling for each of the sites. The interactions corresponding to a given value of / 
in (2.1) couple all the sites within each box of size 2 with the same strength (|) . 
The model has a free parameter c for which we shall use the parametrization 

c = 2^"^. (2.2) 

The parameter of the epsilon-expansion will be defined as 

e = 4-D (2.3). 

This choice has been justified in Ref. [6], but different conventions exist in the 
literature. In the actual calculations reported below, we have selected the values 
and 1 for e. 

We are interested in calculating the magnetic susceptibility of the hierarchical 
Ising model. This quantity can be calculated easily if we know the probability 



for the total spin denoted Pn{S). This probabihty is obviously /^-dependent even 
though we shall not write it explicitly. This probability will be calculated recur- 
sively using the RG method without a rescaling of the spins. We first integrate the 
spins inside boxes of size 2 keeping the sum of the spins in each box constant. We 
then include the terms with / = 1 in (2.1) in a new local measure for these sums. 
We repeat this procedure n times and obtain a measure for the total spin which can 

be normalized as probability. Note that a choice of indices in, j'^i+i completely 

specifies a box of size 2 with the subdivision described above. We call the sum 
of the spin inside this box 5'j„,,....,ij^-^. It can take all the even values between —2 
and 2^ Obviously, S'j„,....,i,_^, = 5'j„,....,i,_^,,i + S'j„,....,i,^,,2- With these notations, the 
recursion formula reads 



22 Pl{Sir„....,il + l,l)Pl{Si^,....,ii+^,2) 

'->i^,....,i^_l_-|^^Oi^,....,i^_l_-|^,lH-Oi^,....,i^_l_-|^,2 



[2A) 



The constant Q+i is adjusted in such a way that the sum of the probabilities 
add up to 1. Strictly speaking, it is not necessary to impose such a normalization 
during the intermediate steps of the calculation, however it keeps the numbers 
reasonably small. This recursion formula has been implemented with a computer 
program. In order to calculate Pn{S), we have to repeat 2"" times a calculation 
involving roughly 2" operations. Consequently, the time necessary to calculate 
Pn{S) scales approximately like 4". With the fastest computer at our disposal, a 
DEC alpha 3000/400, it takes about 10 minutes to calculate Pn{S) from Piq{S) 
when programmed in FORTRAN. 

Due to the size of the calculation, it is clearly necessary to check for round-off 
errors. We have studied the size of these errors by repeating the calculation with 
double-precision instead of simple precision for a large sample of values of (3. For 
n up to 16, we have found very good agreement between the two calculations, the 
differences showing up in the sixth significant digit of the susceptibility. Differences 



in the fifth digits were observed for n = 17 and 18. Calculations for n = 19 and 
beyond have shown less stability, require a lot of computer time and will not be 
reported here. 

The main quantity of interest for us is the average of the square of the total 
spin denoted X„(/3) and defined as 

Xn{P)=J2Pn{S)S^ (2.5) 

s 

Above the critical temperature, this quantity divided by the number of sites has 
a finite thermodynamic limit. It will be called the magnetic susceptibility per site 
and denoted Xn{P) 

XniP) = ^ (2.6) 



3. The Scaling Laws and the Determination of f3c 

The RG method has definite predictions for the the large n behavior of X„. 

.—-. [81 

This is explained at length for instance in chapter 7 of Parisi's textbook. In the 
following we shall just recall the main results and show that our numerical calcula- 
tions reproduce these results with good precision. The case D = 3 will be discussed 
in complete detail, while the results for D = 4 will only be briefly commented upon. 

It is convenient to express X^ as a power of the square of the number of sites, 
namely 

Xni(3) = 22"'^(^'") (3.1) 

The value of u!{l3,n) can be easily estimated for large or small values of (3. For 
large values of f3, all the spins tend to align and uj{f3, n) gets close to 1. For small 
values of /3, the spins at different sites become uncorrelated and the variance of the 
total spin can be approximated by the sum of the individual variances. In other 

f8l 

words ij{P,n) get close to 1/2. A more detailed analysis shows that for large n. 



ijj{f3,n) is attracted by 1 if /5 exceeds a critical value denoted j3c, and by 1/2 if j3 
is less than f3c- When (3 is exactly Pc, oo{(3, n) tends to ^ + ;d • This can be derived 
from the fact that if we reabsorb a factor a/c/4 = 22 + d in the spin variable after 
the integration described in the previous section, the Hamiltonian is invariant in 
the thermodynamic limit. 

We shall now illustrate these results and calculate /3c for D = ?>. In order to get 
a rough idea of the value of /9c, we have plotted in Fig. 1 the trajectories uj{f3,n) 
from n = 1 to n = 16 and for f3 = 0.1, 0.2, ..., 1.6. It appears clearly that the 
separation between the two domains of attraction occurs for a value of [3 between 

1.1 and 1.2. We have repeated this calculation for (3 = 1.10,1.11, ,1.20 as 

shown on Fig. 2. This restricts f3c to the interval [1.17, 1.19]. At this point, it is 
more informative to consider ratios of susceptibilities at two successive values of 
n, because this quantity depends less sensitively on the constant of proportionality 
appearing in the scaling laws. A short calculation shows that 

lim„^oo^o^2 7-T^ = 7: 3.3 

Xn{Pc) D 

Fig. 3 displays -^0^2 ( '^"^(i) ) for 3 values of (3 and for n up to 18. For (3 = 1.179, 
this quantity stays within one percent of the critical value 2/3 for the last six 
iterations. If (3 is increased or decreased by 0.001, the same quantity, departs 
by about 10 percent from 2/3 when n = 18. From this, we conclude that 1.179 
approximates (3c with a precision better than 0.001. 

Applying the same procedure for D = 4, we found that jSc should be inside 
the interval [0.66,0.67]. Fig. 4 shows that 0.665 approximates f3c with a precision 
better than 0.001. 



4. Fitting the numerical data with a simple power law 

In this section, we report our numerical results concerning the /5-dependence 
of the magnetic susceptibility per site at -D = 3 and D = 4 and for f3 < f3c- Again, 
we start with the case D = 3 and discuss it in detail while the case D = 4 will be 
presented more rapidly later in this section. 

We have calculated XniP) for D = 3 and n up to 16 and for values of /3 between 
and 1.18 separated by intervals of length 0.01. As expected, the susceptibility 
rises sharply when f3 gets close to 1.18. We have displayed the results for n = 16 
and (3 < 1.1 on Fig. 5. Results for n = 14 or n = 15 would have been hardly 
distinguishable from n = 16 on a graph of this size. On the other hand, the n- 
dependence becomes more sizable when /3 is closer to its critical value as shown on 
Fig. 6 which also includes results for n = 17. Except for the difference of scales, 
the resemblance between Fig. 5 and Fig. 6 is striking. This can be understood 
from the fact that the numerical data can be fitted very precisely with a simple 
power law of the form 

Xn{(3) = {l-(3/(3c)-^ (4.1) 

in the whole interval [0, 1.18]. In order to justify this claim, we first notice that (4.1) 
implies the inverse logarithmic derivative of the susceptibility is a linear function, 
namely 

i^Logixnim-' = {l)-\Pc-f3) (4.2) 

We now approximate this inverse logarithmic derivative by 

A,/ALo„x„(;3.A,/2)) ^ ,,^(^„(,,^^_,,^(^„(,)) (4^3) 

with A/9 = 0.01, the interval used here. This function is plotted in Fig. 7 for 
n = 16. Remarkably, the numerical data is barely distinguishable from the least 



square linear fit 

Ap/ALogixwiP)) =0.80839-0.67843/5 (4.4) 

The graphs for n=14 or 15 are almost identical. The linear fits are respectively 

Ap/ALog{xi5{P)) =0.80849 - 0.67712/? 

(4.5) 
A(3/ALog{xu{(3)) =0.80867-0.67508/3 

In order to see the corrections to the linear behavior, we have plotted in Fig. 8 
the difference between the linear fit and the data, denoted E{f3), for n = 14, 15 
and 16. These difference are not larger than 0.003 in absolute values and have 
interesting regularities. The reasons why the power law (4.1) is so accurate and 
the nature of its corrections are being investigated. 

We have intentionally used simple precision data to plot Fig. 8 in order to 
give an idea of the round-off errors of the method. Small irregularities are visible 
especially in the low /3 region where their typical size is 10^ . On the other hand, 
E{f3) is smoother for f3 > 0.8. The size of these numerical errors is compatible with 
the claim made in section that numerical errors affects only the sixth significant 
digit of the susceptibility. Indeed, if one of the Xn is replaced by Xn(l + 5) in Eq- 
(4.3), this creates an error 5 A(3 / {ALog{x))'^ ■ A simple inspection shows that 5 can 
be amplified by a factor of order 100 if /5 is not too close to /3c. 

A similar procedure has been followed for D = A. In Fig. 9 we desplay the 
quantity Aj3 / ALog{xi&{l3)) ■ Again the departures from linearity are small and 
the changes in the coefficients of the linear fit for n = 14 and 15 are of the same 
order as those for D = ?>. The difference between the linear fit and the data for 
n = 14, 15 and 16 is displayed on Fig. 10. 



5. Comparisons with High- Temperature and Epsilon Expansions 

Approximate values of the critical exponent 7 can be extracted from linear fits 
presented in the last section. For instance, identifying Eqs.(4.2) and (4.4) yields 
7 = 1.474 and /3c = 1.192. Changes of 1 or 2 percent can be obtained if we restrict 
j3 to intervals closer to I3c- These results can be compared with analytical results 
obtained using the high temperature or the e-expansion. 

The high-temperature expansion of the susceptibility reads 

Xn{l3) = 1 + 6(l,n)/3 + &(2,n)/3^ + - (5.1) 

A straightforward but tedious calculation yields 

\i,n) = (1 - ^r\-Ai - (^)")(i - ^7;)-' - (2'^ - i)(7r') (5.2.1) 



4/ v^v V2' '-^ 2' ^ 'M' 



and 



(5.2.2) 



h2,n) =ihi,n)f - (1 - l)-%n^ - ijr)ii - j)-' 

- 2(^)-+2(l _ (^)n)(i _ |)-1 _ (2- _ 1)(^)2(-+1)] 

Comparing Eq.(5.1) and Eq.(4.2) expanded about 0, we obtain 

7 = &^-ir' (5.3) 

On the other hand, the critical exponents of the hierarchical model near Wil- 
son's non-trivial fixed point can be calculated using the e- expansion. With the 
convention of Eq.(2.2), we obtain at first order in e. 

7 = l + 0.1667e (5.4) 

These analytical results are compared with the numerical results in Fig. 11 and 
12 for D = 3 and D = 4 respectively. It appears clearly that the data agree within 

10 



a few percent with the high-temperature calculation performed above, while a 
significant departure from the e-expansion is shown. This is quite clear in the case 
D = 4 where the e expansion predicts 7 = 1. This is less clear in the case D = 3 
where we do not know at which order we need to truncate what is presumably an 
asymptotic series. 

6. Conclusions 

We have studied the magnetic susceptibility per site of the hierarchical Ising 
model. We found that the /9-dependence of this quantity could be fitted quite 
accurately with a simple power law. It would be quite interesting to understand 
the origin of this result. Small corrections to this power law have been observed 
and it might be possible to handle this problem using Callan-Symanzik's equations. 
In any case, it is remarkable that one is able to obtain an accurate information 
about the critical behavior of a model just by calculating two coefficients of a high 
temperature expansion. 

We insist on the fact that we have made no approximations in our numerical 
calculations and that the results presented are exact up to round-off errors. We 
have analyzed these errors and found that they do not affect our conclusions. The 
fact that we were able to reproduce accurately well-understood analytical results 
such as the scaling laws at the critical temperature and the high temperature 
behavior of the susceptibility seem to rule out errors in implementing Eq. (2.4) 
numerically. 

The discrepancy with the e-expansion is quite surprising. Our limited data 
cannot rule out the possibility that this discrepancy decreases very slowly with 
n. Another possibility is that by starting with Ising spins and applying the RG 
transformation, we never get close enough to the region (in the space of theories) 
where the e-expansion is legitimate. This second possibility seems supported by 
preliminary results obtained in the field theory approach of this problem 
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FIGURE CAPTIONS 

Fig. 1 : uj{n, (3) versus n for D = 3 and values of {3 going from 0.1 to 1.6 by steps 
of 0.1 

Fig. 2 : oj{n^ 13) versus n for D = 3 and values of (3 going from 1.1 to 1.2 by steps 
of 0.01 

Fig. 3 : Log2 of the ratio of two successive values of the susceptibility for D = 3 

and for /3= 1.178, 1.179 and 1.180. 

Fig. 4 : Log2 of the ratio of two successive values of the susceptibility for Z^ = 4 
and for /3= 0.6645, 0.6655 and 0.6665 . 

Fig. 5 : The magnetic susceptibility for D = 3, n = 16 and values of (3 going from 
to 1.18 by steps of 0.01. 

Fig. 6 : The magnetic susceptibility for D = 3, n =14, 15, 16 and 17, and values 
of (3 going from 1.00 to 1.18 by steps of 0.01. 

Fig. 7 : A discrete version of the inverse logarithmic derivative of the susceptibility 
and its linear fit for D = 3 and n = 16. 

Fig. 8 : Difference between the discrete version of the inverse logarithmic derivative 
of the susceptibility and its linear fit for D = 3 and n = 14, 15 and 16. 

Fig. 9 : A discrete version of the inverse logarithmic derivative of the susceptibility 
and its linear fit for D = 4 and n = 16. 

Fig. 10 : Difference between the discrete version of the inverse logarithmic deriva- 
tive of the susceptibility and its linear fit ior D = 4 and n = 14, 15 and 16. 

Fig. 11 : Comparison between the high temperature calculation and our numerical 
estimation of the critical exponent 7 for D = 3 and n between 5 and 16. The solid 
line is the result for the e-expansion at first order. 

Fig. 12 : Comparison between the high temperature calculation and our numerical 
estimation of the critical exponent 7 for D = 4 and n between 5 and 16. The solid 
line is the result for the e-expansion. 
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